Transport phenomena in the asymmetric quantum multibaker map 
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By studying a modified (unbiased) quantum multibaker map, we were able to obtain a finite 
asymptotic quantum current without a classical analogue. This result suggests a general method for 
the design of purely quantum ratchets, and sheds light on the investigation of the mechanisms leading 
to net transport generation by breaking symmetries of quantum systems. Moreover, we propose the 
multibaker map as a resource to study directed transport phenomena in chaotic systems without 
bias. In fact, this is a paradigmatic model in classical and quantum chaos, but also in statistical 
mechanics. 
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I. INTRODUCTION 

In recent years many works have investigated different 
kinds of transport phenomena in periodic dynamical sys- 
tems having no external net force or bias (the so-called 
ratchet effect) [1]. This interest is of fundamental char- 
acter, but it is also motivated by the fact that many pos- 
sible applications exist. For example, they can be use- 
ful to understand and develop rectifiers, pumps, particle 
separation devices, molecular switches, and transistors. 
Also, there is great interest in biology, since the work- 
ing principles of molecular motors can be explained in 
terms of ratchet mechanisms [2]. Finally, we would like 
to mention cold atoms and Bose-Einstein condensates as 
promising fields of application, thanks to recent devel- 
opments of the techniques needed to manipulate them 
[3, 4]. In the following we will have in mind systems that 
present chaotic features, since ratchets generally behave 
this way [5, 6, 7]. Hence, methods from classical and 
quantum chaos become extremely useful. 

The explanation of the appearance of a net current, 
i.e., average momentum different from zero, is one of the 
main topics of this research. In a classical context, this 
behavior can be understood in terms of broken symme- 
tries. One of the most convincing points of view up to 
now, is that all symmetries of the system leading to mo- 
mentum inversion (i.e., sign change) should be broken in 
order to have a net directed current [7] . This amounts to 
say that, being only a necessary condition, we can follow 
Curie's principle, and assume that if the current is not 
forbidden by symmetries, then it should be present. In 
non-Hamiltonian cases, the asymmetrization of a chaotic 
attractor leads to a net directed current [5]. In Hamil- 
tonian systems the asymmetrization of a chaotic layer 
embedded in a mixed phase space, for example, has the 
same consequences for a set of initial conditions inside of 
it [8, 9]. In this case, a mixed dynamics is in general a 
necessary condition (for a notable exception see [10, 11]). 

The vast majority of the previously mentioned papers 
were focused on the classical aspects, leaving the quan- 
tum side less explored. However, there has been sev- 



eral recent publications that deal with this part of the 
problem. These works regard both the experimental [12] 
(systems of cold atoms, mainly) and the theoretical sides 
[8, 10, 13]. In general the quantum versions share sym- 
metry aspects with their classical counterparts, showing 
the corresponding current. But sometimes the relation 
between symmetries and the generated current are less 
obvious in the quantum case. Tunneling, for example, 
can modify the direction of the current [14]. Interfer- 
ence, in any of its forms, generates more complex be- 
haviors [15]. In fact, we will see that the quantum and 
classical behaviors can be very different. 



A remarkable situation appears in some cases, when 
one can find a net quantum current that does not have 
a classical counterpart. This was essentially studied in 
Hamiltonian (non dissipative) systems. The first time 
this phenomenon was found was in the modified kicked 
rotor (KR) at quantum resonance (i.e., the usual KR with 
a second harmonic in the kick) [16, 17]. This is a clas- 
sically chaotic system where the time reversal symmetry 
changing the sign of the momentum is kept. Later, in 
the case of the modified kicked Harper model [18], the 
classical current was found to be exceedingly small in 
comparison with the quantum one. In all of these cases 
the current does not reach an asymptotic value, and in 
fact these systems were called quantum ratchet accelera- 
tors. We present here a different behavior by means of a 
modified multibaker map [19, 20]. This system models a 
particle evolving with free flights and collisions with fixed 
scatterers. We were able to find a finite asymptotic cur- 
rent that is only present in the quantized version. This 
suggests a general method for obtaining purely quantum 
ratchets (no classical current), without unbounded ac- 
celeration. We propose this system as a model to study 
general quantum transport phenomena in the presence of 
asymmetries. 
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II. THE MODEL 



term is defined as 



The well known classical baker's transformation is an 
area-preserving map defined on the unit square phase- 
space (0 < < 1). Here we use an asymmetric version 
which divides the phase space into two regions with dif- 
ferent areas and Lyapunov exponents. The general form 
of this map can be written in terms of one parameter 

s e (0, 1) 



Bs{q,p) 



{lq,sp) 0<q<s 
^{il-s)-\q-s),{l-s)p + s) s<q<l 

(1) 

where we can recover the usual symmetrical case by set- 
ting s = 1/2. The geometric action of the map is repre- 
sented in Fig. 1 (top). 
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FIG. 1: Geometric action of the asymmetric baker's map on 
each cell (top) and of the composition with the translation, 
i.e., the asymmetric multibaker map (bottom). 

The quantum version of the map is defined in a discrete 
/^-dimensional Hilbert space with h = 1/D in terms of 
the quantum Fourier transform with antiperiodic bound- 
ary conditions [21, 22]. 
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where the allowed values of s are such that Di = sD and 
D2 = D — Di are integer numbers. The same procedure 
can be done to obtain an entire family of asymmetric 
quantum baker maps (QBMs) [23, 24]. 

The classical multibaker map [19] is defined in a two- 
dimensional lattice where the phase space at each lattice 
site is a unit square. The dynamics of the map is a com- 
bination of transport to neighboring cells and a local evo- 
lution within a cell. The map is defined as Ms = BgoT, 
where the baker term is the asymmetric baker map de- 
fined in Eq. 1 applied on each cell m, and the transport 
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(m 
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< g < 1/2 
1/2 < g < 1 



(4) 



The geometric action of the multibaker map in the phase 
space of a lattice of squares is represented in Fig. 1 (bot- 
tom). While the baker's map is asymmetric, the trans- 
port term is unbiased as phase space volume is trans- 
ported symmetrically, as can be clearly seen. Notice that 
the transport is entirely due to translations and therefore 
there are no tunneling effects from cell to cell. 

This transformation is not biased neither in the p nor 
in the q coordinate; the reason is that the baker trans- 
formation maps the unit square onto itself and the trans- 
port term is balanced, as previously explained. This is 
the equivalent to the zero average net force typical of 
the dynamical systems addressed in studies of directed 
transport [1]. In what follows we will focus on the coarse 
grained current, that in the classical case is defined by 



Jciassit) = {m{t)) - {m{t - 1)) 



(5) 



In this expression {m{t)) is the average value of the cell 
position m for a given ensemble of initial conditions, at 
a time t. It is easy to see that this definition does not 
take care of the fluctuations inside of each cell. However, 
we underline that there is no bias, making this model 
completely general. 

From the classical point of view, the presence or ab- 
sence of an asymptotic current follows the general criteria 
specified in [7] and depends on the symmetries that re- 
verse the sign of the transport. Here there are two such 
symmetries 



Si: 
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(6) 
(7) 



The first one maps Eg to and therefore is broken 

unless s = 1/2. Sn is a time reversal symmetry and is 
preserved at all times. Thus, according to the criteria of 
[7] there cannot be an asymptotic classical current for un- 
biased initial conditions. Transient effects can be present 
for biased conditions but will die off very rapidly due to 
the exponential mixing property of the Baker map. Nu- 
merical calculations confirm this expectation (see inset 
of Fig. 2(b) and caption). 

The quantum multibaker map is the composition of a 
translation on the lattice site that depends on the value 
of the projectors acting on the right and left subspaces 
of the baker map at each cell [20]. This can be written 
as 



M, 



B.^oT 



^Pl). (8) 



where Pr and Pl are the projectors and t/ is a uni- 
tary translation operator acting on the lattice subspace 
U\m) = |m + 1) (with {|m), m = . . . , -2, -1, 0, 1,2,.. .} 
the basis set of the lattice). 
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III. NUMERICAL RESULTS 

The discrete time propagation of an initial state po 
Plat <^ Pqbm is 



(9) 



Here we will focus on initial states localized in one site 
of the lattice piat = |0)(0|. A word regarding the anal- 
ogy with quantum walks is in order here. In fact, this 
system can be thought as the D-dimensional quantum 
baker map coupled to a quantum walker in an infinite di- 
mensional lattice [25, 26]. The role of the quantum coin 
in the quantum walk is performed here by the quantum 
baker map. The fact that the coin is classically ergodic 
is however an important difference, and is a fundamental 
reason for having no classical current. 

The transport properties of the system can be com- 
puted with the coarse grained density of probabilities of 
the lattice. This distribution is obtained by tracing out 
the internal degrees of freedom inside each cell and pro- 
jecting on the lattice basis. Being 

P(m, t) = Tr [(|m) (m| /qbm) p{t)] (10) 

the mean value of the coarse grained position and the 
quantum current become 

CO 

{x{t)) = Yl mP{m,t) (11) 

171= — CO 

J{t) = {x{t))-{x{t-l)). (12) 

We now turn to analyze the current behavior of this 
system by means of numerical simulations. First, we 
point out that if we take a mixed superposition of cen- 
tered eigenstates of momentum as the initial condition 
for the QBM, we find that there is no current (J) for 
the symmetric case {s = 1/2). This confirms the need to 
break the Si symmetry in order to have a net directed 
transport. Throughout the following calculations we con- 
sider as initial condition the mixed state corresponding 
to the incoherent superposition of the two central mo- 
mentum eigenstates, i.e. 



PQBM = -j^D 



D 
~2 



D 
~2 



G 



(13) 

Fig. 2(a) shows the coarse grained current as a func- 
tion of time with Hilbert space dimension D = 32. 
We have plotted the results corresponding to its maxi- 
mum negative and positive asymptotic values, for which 
Di = 17 and Di = 30 respectively {s = Di/D G [0.5, 1)). 
It is worth mentioning that J rapidly reaches a station- 
ary behavior (at about t = 50) with fluctuations centered 
around a finite non-zero value as it is illustrated by the 
smoothed current in Fig. 2 (solid lines). This shows that 
no unbounded acceleration is present, so there is a recti- 
fication of transport rather than an effective force in this 
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FIG. 2: (Color online) (a): Coarse grained current J as a 
function of time for the QBM with D = 32 and Di = 17 
(dot-dashed line) and Di = 30 (dotted line). In solid lines 
the smoothed current over 20 steps is shown. The solid line at 
J = corresponds to the classical current calculated for 10^ 
initial conditions, (b): The probability distribution P(m, t) 
at time t = 200 for the same cases as before, s = 17/32 and 
s = 30/32 in the classical and quantum version with D = 32. 
In the inset we show how the net current become null in the 
classical case for different initial conditions, calculated for 10^ 
points in phase space, in particular for two squares of area 
1/16 centered in {q,p) = (1/8, 1/2) (s = 0.6 (+) and s = 0.8 
(x)) and (3/8, 1/2) (s = 0.6 (*) and s = 0.8 (□)). 



purely quantum ratchet. The classical current for ana- 
log initial conditions is also displayed, being zero at all 
times. In Fig. 2(b) we show the probability distribution 
P(m, t) and its classical version. 

In Fig. 3(a) we can see the averaged (asymptotic) 
coarse grained current (J) as a function of the Hilbert 
space dimension D for Di = D — 2 and Di = D/2 — 1 (in 
order to make this averages we have taken the values of 
the current from t = 100 up to t = 450). These cases are 
those at which J approximately reaches its maxima in ab- 
solute value (positive and negative current, respectively). 
We observe a smooth dependence of this quantity on D 
(apart from small fluctuations), showing that the effect 
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is generic. In fact, this situation is completely different 
from previous cases [16, 17]. In these works a purely 
quantum current was found only at resonant values of a 
modified KR. 

In Fig. 3(b) we show the averaged current vs s for dif- 
ferent fixed values of D with s > 0.5. We do not show 
values for s < 0.5 since J turns out to be an odd func- 
tion of s around s = 0.5 {{Js) = — (^i-s))- In fact, if 
we apply the symmetry transformation Sj to Eq. 9, and 
then trace out the internal degrees of freedom inside of 
each cell we obtain that Ps{m,t) = Pi-s{—m^t) for all 
t. This result is valid for any initial pqbm symmetrical 
under Sj. Therefore, the symmetry of the current be- 
comes clear (details will be presented elsewhere [28]). As 
a consequence, there is a simple way to obtain current 
inversion in our system. 




IV. CONCLUSIONS 



In summary, we have found a finite asymptotic current 
in a modified and unbiased quantum multibaker map; 
this current has no classical counterpart. Moreover, we 
can invert its direction by exploiting a symmetry with 
respect to the Sj transformation. This is much in the 
same spirit as the corresponding current inversions found 
in the literature [6, 11]. 

The behavior of our system is rather similar to what 
happens in quantum walks. Different combinations of 
the initial states and/or the unitary operators associated 
to the quantum coin produce a different bias in the final 
distribution of probabilities [27]. We believe that this 
property, that is of pure quantum origin is the reason for 
finding a net transport in our system. In fact, we think 
that the existence of a net current is due exclusively to 
interference effects which persists in the D ^ oo limit. 
To obtain the classical behavior, i.e., a vanishing current, 
some noise has to be introduced, just as in the case of 
quantum walk [26] . 

We underline that this is a completely different case 
than that explained in [13], where the quantum current 
appearance is explained by means of the desymmetriza- 
tion of Floquet states, fact that can be directly related 
to the corresponding classical properties of the system. 
A detailed study of this will be presented elsewhere [28]. 
Finally, we consider that the features of this map can 
be exploited to design generic Hamiltonian systems that 
behave in a similar way. In that case, there will be the 
possibility to implement them in cold atoms experiments. 



FIG. 3: (Color online) Averaged coarse grained current (J): 
(a) for different values of the dimension of the Hilbert space 
(D) of the QBM with fixed Di = D /2 + 1 smd Di = D - 2; 
and (b) for the rational values of s = Di/D G [0.5; 1) with 
D = 16, 32, 64. The average of J is computed up to 450 
iterations of the different maps beginning with t — 100 (where 
the current approximately reaches its asymptotic behavior). 
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